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ABSTRACT 

Pressure  and  velocity  data  were  collected  in  a  full  scale  jet 
engine  test  cell  in  order  to  validate  the  predictive  accuracy  of  a 
two-dimensional  and  axisymmetric  primitive  variable  computer  code.  It 
was  found  that  the  model  reasonably  predicted  the  velocity  profiles  in 
the  augmentor  tube.  Inaccuracy  increased  at  higher  engine  thrust  settings 
at  positions  far  downstream  in  the  augmentor  tube.  Predicted  pressure 
profiles  were  reasonable  but  the  magnitudes  were  in  considerable  error 
at  high  flow  rates. 
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I.  INTRODUCTION 


Jet  engine  test  cells  are  in  wide  use  today  in  military  and  com¬ 
mercial  aviation  maintenance  programs.  Their  purpose  is  to  provide 
for  testing  of  jet  engines  under  conditions  that  approximate  the  in¬ 
stalled  environment,  thus  enabling  the  responsible  authority  to  veri¬ 
fy  the  engine's  capabilities  throughout  its  operational  envelope 
prior  to  returning  it  to  service.  However,  due  to  exhaust  pollution 
problems  inherent  to  their  operation,  the  use  of  these  test  facilities 
has  recently  encountered  a  considerable  amount  of  criticism  and  even 
judicial  action  in  the  State  of  California. 

The  necessity  for  testing  an  engine  after  overhaul  is  apparent. 
However,  the  pollution  problems  that  result  when  testing  today's  high 
power,  high  mass  flow  engines  in  present  test  facilities  must  also  be 
considered. 

"The  typical  test  cell  incorporates  an  inlet,  a  horizontal  test 
section  and  vertical  exhaust  stack.  The  engine  to  be  tested  is 
normally  mounted  near  the  center  of  the  cell  to  allow  the  development 
of  a  nearly  uniform  engine  inlet  velocity  profile.  The  engine  ex¬ 
hausts  into  an  augmentor  tube  which  entrains  additional  air  for  ex¬ 
haust  gas  cooling  and  dilution.  The  quantity  of  this  secondary  air 
is  crucial  to  proper  engine  testing  and  test  cell  performance" 

[Ref.  1].  The  ratio  of  entrained  air  mass  flow  rate  to  engine  mass 
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flow  rate  is  known  as  "augmentation  ratio".  An  augmentation  ratio  that 
is  too  high  could  result  in  inaccurate  engine  performance  measurements 
and  in  exceeding  the  test  cell  structural  limits  due  to  pressure 
gradients  and  excessive  cell  pressure  reduction.  Insufficient  secondary 
air  could  allow  engine  exhaust  gas  recirculation  to  the  engine  inlet  and 
cause  excessive  opacity  of  visible  emittants.  It  is  therefore  necessary 
to  develop  mathematical  models  which  can  be  used  to  predict  the  flow 
fields  within  turbojet  test  cells.  These  models  could  then  be  used  as 
a  cost  effective  means  for  determination  of  optimum  test  cell  designs 
for  meeting  the  military  engine  testing  needs  while  minimizing  noise  and 
chemical  pollution  effects  on  the  local  environment. 

With  this  in  mind,  considerable  experimental  and  computer  modelling 
work  has  been  done  at  the  Naval  Postgraduate  School  to  further  the 
understanding  of  the  flow  field  and  detailed  operating  characteristics 
encountered  in  typical  test  cells.  This  knowledge  is  currently  needed 
to  facilitate  the  necessary  modifications  to  present  test  cells  and  for 
the  design  of  new  cells  to  accommodate  the  high  technology  engines  of 
the  future. 

The  two-dimensional  modelling  of  the  flow  within  turbojet  test 
cells  was  initiated  at  NPS  by  Hayes  and  Netzer  [Ref.  2].  They  adapted 
the  two-dimensional,  ip-oj  model  of  Spalding,  et  al  [Refs.  3,  4]  to  the 
test  cell  geometry  and  operating  conditions.  This  work  was  continued 
by  Speakman,  Hayes  and  Netzer  [Ref.  5].  The  model  made  use  of  the 
fundamental  equations  of  conservation  of  mass,  momentum  and  energy. 

The  dependent  variables  were  stream  functions  (tp) ,  vorticity  (co). 
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turbulence  kinetic  energy  (K),  the  turbulence  dissipation  rate  (e),  and 
the  temperature  (T).  As  a  consequence  of  the  ip-oj  formation,  pressure 
(and  to  a  large  extent  velocity)  is  removed  from  the  governing  equa¬ 
tions.  These  second  order  elliptic  partial  differential  equations 
were  reduced  to  "finite-difference,"  non-linear  algebraic  equations  and 
solved  simultaneously  with  an  iterative  procedure.  Results  of  this 
computer  simulation  appeared  to  be  quite  accurate  at  low  subsonic  en¬ 
gine  exhaust  Mach  numbers  (M  <  0.6),  but  pressure  predictions  were  un¬ 
reliable  at  high  subsonic  or  sonic  exhaust  velocities.  These  compari¬ 
sons  between  predictions  and  experiment  were  made  by  Walters  and  Netzer 
[Ref.  6]  using  a  one  eighth  scale  model  of  an  NAS  Alameda  test  cell. 
Other  problems  related  to  the  ip-u  formulation  were:  1)  restriction  to 
constant  density  flows,  or  to  flows  in  which  density  varied  only  with 
temperature;  2)  boundary  conditions  were  difficult  to  specify;  3)  con¬ 
siderable  difficulty  was  experienced  in  obtaining  converged  solutions, 
especially  for  non-uniformly  spaced  grids  and  high  flow  rates,  and  4) 
the  ip-co  model  is  not  easily  extended  to  three  dimensional  flows. 

To  help  alleviate  the  above  modelling  difficulties  Stevenson  and 
Netzer  {Ref.  ^utilized  a  primitive  variable  ( pressure- velocity)  model . 
This  computer  model  (adapted  from  the  CHAMPION  2/E/FIX  computer  program 
developed  by  Pun  and  Spalding  ^Ref.  7^  extended  the  modelling  capabili¬ 
ties  to  more  complex  geometries  and  to  high  subsonic  and  sonic  engine 
exhaust  velocities.  It  also  is  readily  applicable  to  variable  density 
flows.  The  earlier  results  of  the  i/j-uj  model  and  those  of  the  u-v-p 
model  were  compared  to  the  empirical  data  obtained  with  the  subscale 
test  cell  {jtef.  6*].  Good  agreement  between  the  two  methods  was  found 
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for  low  subsonic  engine  exhaust  conditions.  The  primitive  variable 
computer  model  produced  reasonable  results  for  the  flow  field  downstream 
of  the  engine  exhaust  up  to  the  sonic  exhaust  condition. 

The  primary  objective  of  this  investigation  was  to  adapt  Spalding's 
CHAMPION  2/E/FIX  primitive  variable  computer  model  (as  modified  by 
Stevenson)  to  the  geometry  and  flow  conditions  of  an  actual  operational 
full  scale  turbojet  test  cell  and  to  test  that  model  with  empirical  data 
from  the  same  cell.  The  turbojet  test  cell  chosen  was  one  located  at 
NAS  Alameda,  operated  by  NARF  Alameda.  The  engine  being  operated  at 
the  time  of  data  acquisition  was  the  General  Electric  TF34-GE-100,  cur¬ 
rently  used  on  the  USAF  A-10  attack  aircraft.  The  major  computer  model 
modifications  required  were:  1)  to  increase  the  geometry  to  the  full 
scale,  and  2)  to  change  the  uniform  engine  exhaust  velocity  profile  to 
a  non-uniform  profile  which  adequately  represented  the  primary  and 
secondary  flows  from  the  TF-34  turbofan  engine. 
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II.  EXPERIMENTAL  INVESTIGATION 


A.  INTRODUCTION 

The  turbojet  test  cell  chosen  for  this  investigation  was  test  cell 

No.  13  located  in  Building  372  at  NAS  Alameda.  It  is  operated  by  NARF 

Alameda  as  the  prime  means  of  evaluating  the  USAF  TF34-GE-100  turbofan 
engine  following  required  normal  periodic  maintenance. 

The  TF34-GE-100  is  a  turbofan  engine  in  current  use  on  the  USAF 
twin  engine  A- 1 0  attack  airplane.  The  engine  is  2.54m  long  and  produces 
a  maximum  rated  thrust  of  9065  lbf  JjRef .  8j.  The  bypass  ratio  is  6.22:1 
with  maximum  rated  airflow  at  151.0  kg/sec.  The  exhaust  nozzle  is  a 
fixed  area,  convergent  type  with  the  cold  stream  being  0.474  sq.  m.  in 
area  and  the  hot  stream  0.148  sq.  m.  in  area.  While  operating  at  maxi¬ 
mum  power  the  core  exhaust  gas  is  estimated  to  have  a  velocity  of  305  m/ 

sec  and  a  temperature  of  393°C  jjRef.  8^. 

To  provide  the  data  for  model  validation  it  was  decided  to  measure 
velocity  profiles  within  the  test  cell  and  at  two  locations  within  the 
augmentor  tube  as  well  as  the  axial  pressure  distribution  throughout 
the  test  cell.  To  most  accurately  determine  velocity  profiles,  mul¬ 
tiple  pitot-static  tubes  with  associated  stagnation  temperature  measure¬ 
ments  would  have  been  required.  However,  due  to  cost  and  time  restraints 
involved  with  the  manufacturing  and  installation  of  such  equipment,  a 
compromise  approach  to  this  problem  was  taken.  Two  stagnation  tube 
rakes,  each  fitted  with  two  thermocouples,  were  installed  in  the  aug¬ 
mentor  tube.  Also,  thirty  two  static  pressure  ports  were  installed  in 
the  augmentor  tube  wall.  Velocity  measurements  were  also  made  within 
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the  test  cell  using  hot  wire  velocimeters .  Details  of  the  data  acquisi¬ 
tion  as  well  as  the  method  of  data  reduction  are  presented  below. 

B.  TEST  CELL 

The  turbojet  test  cell  utilized  in  this  investigation  is  typical  of 
test  cells  currently  used  in  military  and  commercial  aviation  maintenance 
programs.  It  was  a  concrete  blockhouse  type  structure  43m  long,  9m  high, 
and  9m  wide  (Fig.  1).  It  had  a  horizontal  inlet,  which  was  retangular  in 
cross  section,  and  exhausted  into  a  vertical  stack.  Not  shown  in  Figure 
1  are  the  noise  suppression  devices  which  were  located  over  the  cross 
sections  of  both  the  test  cell  inlet  and  back  wall.  Located  15m  from 
the  cell  inlet  was  the  test  platform  which  stood  3m  above  the  floor  of 
the  cell.  The  jet  engine  test  stand  and  the  first  section  of  the  aug- 
mentor  tube  were  located  on  this  platform.  The  TF-34  was  mounted  on  the 
test  stand  with  an  attached,  fixed  geometry  bellmouth  type  inlet,  along 
with  necessary  start-up  and  test  equipment  to  monitor  engine  performance 
in  the  cell  control  room.  This  test  stand  was  movable,  so  that  the 
distance  between  the  exit  plane  of  the  engine  and  entrance  plane  of  the 
augmentor  tube  could  be  varied.  During  this  investigation  this  distance 
was  fixed  at  1 .7m. 

The  first  section  of  the  augmentor  tube  was  located  on  the  test 
platform.  As  shown  in  Figure  1  it  too  was  movable.  It  was  a  steel 
pipe  with  9.5mm  thick  walls  and  was  approximately  1.8m  in  diameter  and 
1.6m  long.  The  next  section  was  also  a  constant  diameter  pipe,  5.2m 
long.  During  this  investigation,  the  pipe  protruded  only  5m  beyond 
the  first  section.  The  radial  gap  between  these  two  sections  was,  on 
the  average,  16mm.  The  next  section  was  3.5m  long  and  diverged  in 
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diameter  from  approximately  1.8m  to  2.7m.  At  the  downstream  end  of  the 
latter  section  was  a  water-injection  ring.  Following  a  0.31m  gap  was 
the  next  section:  a  constant  3.05m  diameter  pipe,  13.1m  long.  In 
this  section  were  three  staggered  rows  of  four  noise  suppression  de¬ 
vices.  Each  unit  was  a  3.2m  long  cylinder,  0.6m  in  diameter. 

C.  INSTRUMENTATION  AND  DATA  ACQUISITION 

Two  pitot  tube  rakes  were  utilized  in  this  investigation.  Each 
rake  was  constructed  from  102mm  wide,  9.5mm  thick,  302  stainless  steel. 
The  cross-like  configuration  used  is  depicted  in  Figs.  2  and  3.  The 
first  rake  was  located  127mm  from  the  entrance  plane  of  the  augmentor 
tube.  The  tips  of  the  pitot  tubes  being  25.4mm  from  the  entrance. 

The  vertical  member  of  this  rake  was  1.67m  high  and  had  seven  pitot 
tubes  built  into  it.  Each  horizontal  member  was  0.822m  long  with  three 
pitot  tubes  inserted.  Horizontal  and  vertical  members  were  attached  as 
shown  in  Fig.  4.  To  hold  the  rakes  in  position,  50.8mm  high  tabs  were 
welded  into  the  augmentor  tube  at  the  12,  3,  6,  and  9  o'clock  positions. 
The  rakes  were  then  attached  to  the  tabs  as  shown  in  Fig.  5.  The  down¬ 
stream  rake  was  constructed  exactly  the  same  as  the  first  except  that 
the  vertical  length  was  1.70m  and  each  horizontal  member  was  0.833m 
long.  This  second  rake  was  installed  5.05m  downstream  from  the  aug¬ 
mentor  entrance  plane  (127mm  upstream  from  the  start  of  the  diverging 
section  of  the  tube).  The  6.35mm  diameter  stainless  steel  tubing  used 
for  the  pitot  tubes  had  a  wall  thickness  of  1.57mm.  Each  tube  pro¬ 
truded  101.6mm  forward  from  the  rake.  The  tubes  were  silver-soldered 
on  the  back  side  of  the  rake  and  routed  from  the  augmentor  tube  at  the 
aft  end  of  each  attachment  tab  through  holes  drilled  in  the  augmentor. 
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Tygon  tubing  (3/16  in.  ID  x  1/16  in.  wall)  was  then  routed  from  pitot 
tube  to  a  Scanivalve  located  in  a  standing  cabinet  next  to  the  jet- 
engine  test  stand.  The  output  of  the  Scanivalve  was  routed  through  the 
floor  of  the  cabinet  and  underneath  the  test  platform  to  the  recording 
apparatus  in  the  control  room  (Fig.  6).  In  addition  to  the  pitot  tubes, 
thermocouples  were  mounted  on  each  rake  as  shown  in  Fig.  2.  The  thermo¬ 
couple  wire  was  routed  from  the  augmentor  with  the  pitot  tubing,  then  to 
ice  bath  references.  The  thermocouple  data  were  obtained  from  a  digital 
voltmeter. 

Static  pressure  ports  were  located  at  the  2,  5,  7,  and  10  o'clock 
positions  at  axial  locations  0.152m,  0.305m,  0.610m,  1.212m,  3.96m, 
5.18m,  6.10m,  and  7.32m  downstream  of  the  augmentor  inlet.  Circumferen¬ 
tial  ports  were  used  to  determine  whether  or  not  the  flow  was  axi sym¬ 
metric.  Tygon  tubing  was  attached  to  Swagelock  fittings  at  each  port 
and  routed  to  a  second  Scanivalve  in  the  standing  cabinet. 

In  order  to  determine  the  velocity  of  the  secondary  (test  cell)  air, 
velocimeters  were  located  0.6m  and  1.8m  above  the  TF34  engine.  The 
velocimeters  used  were  Datametrics  Air  Flow  Meters,  commonly  used  to 
measure  the  velocity  in  air  conditioning  ducts.  Prior  to  their  use, 
calibration  was  performed  at  the  NPS  using  a  small  wind  tunnel.  The 
resulting-  calibration  curve  is  shown  in  Fig.  7.  The  output  of  these 
velocimeters  were  routed  underneath  the  jet  engine  test  stand  to  the 
standing  cabinet  and  then  to  the  control  room  as  described  above.  In 
the  control  room  the  output  signals  were  connected  to  a  digital  volt¬ 
meter. 

The  instrumentation  apparatus  set  up  in  the  control  room  (Fig.  6) 
consisted  of  a  locally  manufactured  Scanivalve  controller,  two  Doric 
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digial  voltmeters  and  two  Honeywell  "Electronic  196",  2  pen  strip  chart 
recorders . 

D.  EXPERIMENTAL  PROCEDURE 

Prior  to  actual  data  acquisition  the  following  steps  were  performed: 

1)  The  30.5  cm  radial  gap  between  the  second  and  third  augmentor 
sections  was  blocked  off  with  1.9  cm  thick  plywood.  This  was 
done  to  prevent  any  secondary  cell  air  from  entering  the  aug¬ 
mentor  downstream  of  the  data  acquisition  equipment.  It  should 
be  noted  that  an  airtight  seal  was  not  attempted. 

2)  All  static  and  dynamic  pressure  ports  were  pressure  checked  to 
insure  proper  connections  and  identification. 

3)  All  thermocouples  were  checked  for  proper  rake  attachment, 
wiring  and  identification. 

The  actual  data  acquisition  process  began  before  anticipated  due  to 
inadequate  thermocouple  installation.  After  initial  TF34  start  up  and 
a  brief  run  time  (during  which  the  recording  instruments  were  calibrated) 
the  test  cell  operators  discovered  an  oil  leak  in  the  TF-34  and  conse¬ 
quently  shut  it  down  for  repair.  This  provided  an  opportunity  to 
inspect  the  test  equipment.  At  this  time  it  was  discovered  that  the 
thermocouples  were  not  suitably  constructed  to  withstand  a  complete 
series  of  tests.  Therefore,  it  was  decided  to  record  termperature  data 
during  the  initial  TF-34  bearing  breakin  tests  at  various  power  settings 
for  as  long  as  the  thermocouples  remained  in  place.  Reliable  tempera¬ 
ture  data  were  only  obtained  for  maximum,  90%,  and  idle  power  settings. 
The  data  appears  in  Table  I. 
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POWER 

SETTING 


100% 


90% 


IDLE 


TABLE  I: 


CORE  FAN  THRUST  RAKE 
RPM  RPM  LBS 

FORWARD 

17228  6667  8611 

AFT 

FORWARD 

16120  5803  5967 

AFT 

FORWARD 

11350  1685  382 

AFT 


THERMOCOUPLE 

TEMPERATURE 

Cc) 

CENTER 

438 

LOWER 

39 

CENTER 

127 

LOWER 

76 

CENTER 

461 

LOWER 

32 

CENTER 

156 

LOWER 

62 

CENTER 

498 

LOWER 

18 

CENTER 

105 

LOWER 

42 

Thermocouple  Data  Obtained  During  Initial  Engine  Testing 
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After  the  brief  bearing  breakin  tests  the  TF-34  was  put  through  a 
series  of  five  minute  duration  tests  at  a  variety  of  power  settings. 
When  a  particular  power  setting  was  reached  two  minutes  were  allowed 
to  insure  steady-state  operation  and  then  the  following  data  were  re¬ 
corded: 

1)  From  the  test  cell  instrumentation: 

a)  ambient  temperature 

b)  ambient  pressure 

c)  cell  pressure 

d)  TF-34  core  RPM 

e)  TF-34  fan  RPM 

f)  TF-34  thrust 

2)  From  NPS  installed  equipment  and  instrumentation: 

a)  static  pressures  along  the  augmentor  wall 

b)  pitot  tube  pressures  at  each  rake 

c)  temperatures  at  each  rake 

Each  48  channel  Scanivalve  was  set  up  so  that  port  number  one  was 
open  to  the  cell  pressure.  A  fixed  pressure  (set  using  a  mercury 
manometer)  was  input  to  port  number  two  in  order  to  provide  a  refer¬ 
ence  calibration  each  time  a  cycle  was  begun.  During  each  run  one 
Scanivalve  was  cycled  using  the  Scanivalve  controller,  data  recorded 
on  the  strip  chart,  connections  switched,  and  the  same  procedure 
followed  for  the  remaining  Scanivalve.  The  pitot  rakes  were  labeled 
as  shown  in  Fig.  8,  and  static  ports  as  in  Fig.  9.  The  recorded  data 
are  presented  in  Tables  II  through  IV,  the  letter  "A"  referring  to  the 
forward  rake,  "B"  to  the  aft  rake,  and  letters  "C  through  J"  to  the 
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static  port  stations  from  front  to  rear.  The  data  obtained  from  the 
intermediate  power  setting  run  was  not  utilized  due  to  lack  of  reliable 
temperature  information  at  that  thrust  level. 

There  were  three  specific  problems  encountered  during  the  data  ac¬ 
quisition  process: 

1)  Failure  of  thermocouple  attachments.  This  resulted  in  lack  of 
temperature  data  during  each  run  that  pressure  data  were  ob¬ 
tained.  However,  this  was  overcome  by  correlating  (by  thrust 
setting  and  fan  and  core  speeds)  the  temperature  data  recorded 
at  the  start  of  the  TF-34  test  series  with  the  pressure  data 
obtained  later. 

2)  During  the  entire  test  series  a  loud  resonant  frequency  was 
heard.  This  initially  forced  engine  shut  down  to  allow  for 
cell  inspection.  No  damage  was  discovered  and  all  instrumen¬ 
tation  appeared  normal  so  testing  was  continued  as  planned. 
However,  upon  removal  of  test  equipment  it  was  discovered  that 
two  of  the  welds  holding  the  aft  rake  in  place  had  broken.  It 
should  also  be  noted  that  at  the  same  time  it  was  discovered 
that  two  pitot  tubes  had  detached  from  the  same  rake.  The 
latter  apparently  occurred  during  the  last  test  since  no 
unusual  change  occurred  in  the  recorded  data. 

3)  Due  to  electrical  connection  difficulties  the  velocimeters  were 
not  installed  during  the  engine  testing  periods  when  other  data 
were  obtained.  When  this  problem  was  corrected  (after  the  pres¬ 
sure  data  were  obtained)  it  was  discovered  that  the  wooden  blocks 
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RUN  NO.:  1 


POWER  SETTING:  100% 


TIME:  1305 


PITOT  RAKE  (P-PCELL)  STATIC  PORT  (P-PCELL) 


SCANIVALVE 

CHANNEL 

PORT  I.D. 

(N/m2xl0-3) 

SCANIVALVE 

CHANNEL 

PORT  I.D. 

(N/m2xl0-3) 

1 

Zero 

1 

Zero 

2 

REF  PRES 

2 

REF  PRES 

3 

A1 

34.30 

3 

Cl 

-14.196 

4 

A2 

39.93 

4 

C2 

-13.776 

5 

A3 

3.56 

5 

C3 

-13.803 

6 

A4 

-0.63 

6 

C4 

-14.266 

7 

A5 

41.11 

7 

D1 

-13.218 

8 

A6 

5.97 

8 

D2 

-12.797 

9 

A7 

-0.48 

9 

D3 

-12.956 

10 

A8 

39.58 

10 

D4 

-13.190 

11 

A9 

3.78 

11 

El 

-9.694 

12 

A10 

-0.84 

12 

E2 

-9.370 

13 

All 

40.19 

13 

E3 

-9.432 

14 

A12 

2.68 

14 

E4 

-9.860 

15 

A13 

-0.48 

15 

FI 

-8.226 

16 

B1 

26.84 

16 

F2 

-8.226 

17 

B2 

19.81 

17 

F3 

-8.226 

18 

B3 

8.39 

18 

F4 

-8.226 

19 

B4 

3.30 

19 

G1 

-6.433 

20 

B5 

23.24 

20 

G2 

-6.661 

21 

B6 

10.36 

21 

G3 

-6.530 

22 

B7 

3.12 

22 

G4 

-6.433 

23 

B8 

18.71 

23 

HI 

-6.268 

24 

B9 

7.51 

24 

H2 

-6.171 

25 

BIO 

1 .28 

25 

H3 

-5.875 

26 

Bll 

15.59 

26 

H4 

-6.33 

27 

B12 

5.76 

27 

11 

-4.144 

28 

B13 

0.04 

28 

12 

-4.013 

29 

13 

-3.820 

30 

14 

-4.047 

31 

J1 

-1.731 

32 

J2 

-1 .958 

33 

J3 

-2.027 

34 

J4 

-1 .827 

FAN 

RPM:  6701 

CORE 

17360 

THRUST:  8685 
(LBS) 

CELL 
TEMP:  17.8 

(°c) 

RAKE  A 

RAKE  B 

CELL  AMBIENT 

PRESSURE:  1 .01 29x10^  1.0141xl05 

(N/m2) 

TABLE  II:  EXPERIMENTAL  DATA,  100%  POWER 
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RUN  NO:  3 

PITOT  RAKE  (P-PCELL 

POWER  : 

) 

SETTING:  90% 

STATIC 

PORT  (p- 

TIME:  1327 

PCELL> 

SCANIVALVE 

CHANNEL 

PORT  I.D. 

(N/m2xlO 

-3x  SCANIVALVE 
>  CHANNEL 

PORT  I.D. 

(N/m2xlO-3) 

1 

Zero 

1 

Zero 

2 

REF  PRES 

2 

REF  PRES 

3 

A1 

22.64 

3 

Cl 

-9.29 

4 

A2 

30.86 

4 

C2 

-10.08 

5 

A3 

0.49 

5 

C3 

-9.56 

6 

A4 

-0.39 

6 

C4 

-9.29 

7 

A5 

29.04 

7 

D1 

-9.23 

8 

A6 

4.18 

8 

D2 

-9.23 

9 

A7 

-0.19 

9 

D3 

-9.33 

10 

A8 

30.60 

10 

D4 

-9.45 

11 

A9 

2.58 

11 

El 

-7.46 

12 

A10 

-0.29 

12 

E2 

-6.70 

13 

All 

31 .71 

13 

E3 

-7.07 

14 

A12 

3.65 

14 

E4 

-7.55 

15 

A13 

-0.26 

15 

FI 

-5.71 

16 

B1 

20.64 

16 

F2 

-5.62 

17 

B2 

16.86 

17 

F3 

-5.53 

18 

B3 

8.18 

18 

F4 

-5.59 

19 

B4 

1.60 

19 

G1 

-4.14 

20 

B5 

17.30 

20 

G2 

-4.41 

21 

B6 

8.41 

21 

G3 

-4.41 

22 

B7 

3.16 

22 

G4 

-4.26 

23 

B8 

13.96 

23 

HI 

-4.08 

24 

B9 

5.83 

24 

H2 

-3.87 

25 

BIO 

1.38 

25 

H3 

-3.75 

26 

Bll 

12.72 

26 

H4 

-4.12 

27 

B12 

3.96 

27 

11 

-2.76 

28 

B13 

0.98 

28 

12 

-2.79 

29 

13 

-2.45 

30 

14 

-2.61 

31 

J1 

-1 .10 

32 

J2 

-1.40 

33 

J3 

-1.40 

34 

J4 

-1.34 

FAN 

CORE 

RPM:  5933 

16403 

THRUST:  6304 

(LBS) 

CELL 

RAKE  A  RAKE  B 

CELL  , 

AMBIENT 

TEMP:  18.3 
(°C) 

PRESSURE:  1 .01 36x1 0D 
(N/nr ) 

1 .0144x10 

TABLE  III:  EXPERIMENTAL  DATA,  90%  POWER 
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RUN  NO:  4 

POWER  SETTING 

:  IDLE 

TIME:  1337 

PITOT 

RAKE  (P“Pq£[[ 

) 

STATIC  PORT  (p- 

PCELL^ 

SCANIVALVE 

CHANNEL 

PORT  I.D. 

(N/m2xl0-3) 

SCANIVALVE 

CHANNEL 

PORT  I.D. 

(N/m2xl0~3) 

1 

Zero 

1 

Zero 

2 

REF  PRES 

2 

REF  PRES 

3 

A1 

1.44 

3 

Cl 

-0.64 

4 

A2 

2.56 

4 

C2 

-0.77 

5 

A3 

0.45 

5 

C3 

-0.74 

6 

A4 

0.27 

6 

C4 

-0.71 

7 

A5 

2.52 

7 

D1 

-0.71 

8 

A6 

0.41 

8 

D2 

-0.74 

9 

A7 

0.27 

9 

D3 

-0.74 

10 

A8 

2.43 

10 

D4 

-0.71 

11 

A9 

0.27 

11 

El 

-0.58 

12 

A10 

0.14 

12 

E2 

-0.49 

13 

All 

2.48 

13 

E3 

-0.52 

14 

A12 

0.23 

14 

E4 

-0.55 

15 

A13 

0.14 

15 

FI 

-0.40 

16 

B1 

1  .53 

16 

F2 

-0.43 

17 

B2 

1 .48 

17 

F3 

-0.40 

18 

B3 

0.54 

18 

F4 

-0.43 

19 

B4 

0.50 

19 

G1 

-0.28 

20 

B5 

0.27 

20 

G2 

-0.28 

21 

B6 

1.35 

21 

G3 

-0.25 

22 

B7 

1.35 

22 

64 

-0.28 

23 

B8 

0.36 

23 

HI 

-0.28 

24 

B9 

1.12 

24 

H2 

-0.18 

25 

BIO 

0.54 

25 

H3 

-0.18 

26 

Bll 

0.32 

26 

H4 

-0.15 

27 

B12 

1.12 

27 

11 

-0.03 

28 

B13 

0.41 

28 

12 

-0.03 

29 

13 

0 

30 

14 

0 

31 

J1 

-0.03 

32 

J2 

0 

33 

J3 

0 

34 

J4 

0 

FAN  CORE 
RPM:  1801  11404 


THRUST:  410 
(LBS) 

CELL  RAKE  A  RAKE  B 
TEMP:  18.3 
(°C) 


CELL  ,  AMBIENT  , 

PRESSURE:  1.0144x10®  1. 0144x10® 

(N/m2) 


TABLE  IV:  EXPERIMENTAL  DATA,  IDLE  POWER 
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at  the  intersection  of  the  second  and  third  augmentor  sections 
had  partially  detached.  This  caused  the  velocity  of  the  cell 
air  to  probably  be  somewhat  higher  during  the  tests  using  the 
velocimeters  than  it  was  during  the  pressure/acquisition  runs 
when  the  blocks  were  in  place.  The  velocimeter  data  are  shown 
in  Table  V.  These  data  were  used  only  to  obtain  the  general 
shape  of  the  velocity  profile  above  the  engine. 


E.  DATA  REDUCTION 

The  (p-pcell)  data  that  appear  in  Tables  II  through  IV  came  directly 
from  the  strip  chart  recordings,  or,  in  the  case  where  the  reading  was 
too  small,  from  the  digital  voltmeter.  The  calibration  scale  on  the 
strip  chart  (psi/division)  was  obtained  from  a  1.05  in.  Hg  reference  in¬ 
put  to  Channel  No.  2  of  both  Scanivalves.  Using  this  information, 
axial  pressure  profiles  along  the  four  circumferential  positions  were 
plotted  and  an  average  pressure  profile  was  then  estimated  for  each 
power  setting  (Figs.  10,  11,  and  12).  These  profiles  were  then  used  to 
obtain  the  static  pressure  at  the  aft  rake  location  (it  was  assumed 
that  pressure  was  uniform  across  the  augmentor  at  the  aft  rake). 

The  measured  stagnation  pressure  at  each  pitot  tube  together  with 
the  static  pressure  at  the  rake  were  used  to  obtain  a  Mach  number  pro¬ 
file  using  the  expression 


X 


Y 
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Velocity  RPM  Thrust 

(m/sec)  (lbs) 


Upper 

Velocimeter 

Lower 

Velocimeter 

Fan 

Core 

0.51 

1.14 

1789 

11265 

520 

2.03 

3.05 

6780 

1  7285 

9072 

TABLE  V:  VELOCIMETER  DATA 
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with  y  assumed  equal  to  1.4.  The  thermocouple  readings  were  plotted  and 
used  to  estimate  an  assumed  axi symmetric  stagnation  temperature  profile 
for  the  aft  rake  at  each  power  setting  (Fig.  13).  The  estimated  stagna¬ 
tion  temperature  profile  and  local  Mach  number  were  then  used  to  deter¬ 
mine  a  static  temperature  for  each  pitot  tube  location.  Using  this  in¬ 
formation  and  the  properties  of  air  at  standard  temperature  and  pressure, 
accoustic  speed  and  velocity  were  found  for  each  pitot  tube  position. 
Static  temperature  for  the  aft  rake  is  plotted  as  a  function  of  radius 
in  Fig.  14,  velocity  profiles  are  shown  in  Figs.  15,  16,  and  17. 

The  static  pressure  of  the  front  rake  was  estimated  by  extrapolating 
the  average  pressure  profile  to  the  augmentor  inlet  position.  This  was 
required  due  to  the  fact  that  the  first  set  of  static  pressure  ports 
was  installed  15.2  cm  aft  of  the  augmentor  entrance  plane.  This  static 
pressure  was  assumed  to  be  uniform  across  the  augmentor  for  each  power 
setting.  Mach  number  profiles  were  then  obtained  as  above.  At  the 
forward  rake  y  was  not  assumed  to  be  uniform  across  the  augmentor.  In¬ 
stead,  y  was  taken  to  be  1.33  at  the  centerline  position  and  at  the  four 
pitot  tubes  located  circumferentially  25.4  cm  from  the  center.  At  the 
remaining  positions  y  was  assumed  to  be  1.4.  This  property  change  was 
made  since  heat  induced  discoloration  on  the  forward  rake  was  observed 
to  occur  from  the  core  jet  to  a  radius  of  approximately  30  cm.  Stagna¬ 
tion  temperature  profiles  were  estimated  using  the  measured  temperatures 
at  the  two  rake  locations  and  are  shown  in  Fig.  18.  Using  this  stagna¬ 
tion  temperature  distribution  with  the  Mach  number  profiles,  accoustic 


27 


speed  and  velocity  were  determined  as  above.  The  resulting  velocity 
profiles  for  each  power  setting  are  shown  in  Figs.  19  through  21. 


28 


III.  COMPUTER  MODEL 


A.  INTRODUCTION 

The  computer  program  used  to  model  the  turbojet  test  cell  of  this 
investigation  was  the  CHAMP  ION/2/ E/F IX  computer  program  developed  by 


This  is  a 


Pun  and  Spalding 


Stevenson 


primitive  variable  (pressure-velocity)  model  for  two-dimensional,  turbu¬ 
lent  recirculating  flows.  The  equations  are  cast  into  finite  difference 
form  and  solved  using  a  line-by-line  iteration  method.  A  detailed  dis¬ 
cussion  of  the  model  can  be  found  in  Refs.  1  and  7. 

B.  PHYSICAL  MODEL 

The  major  modifications  made  to  the  computer  program  as  utilized  by 
Stevenson  ^Ref.  ij  were  to  change  the  geometry  of  the  test  cell  and  to 
incorporate  the  engine  exhaust  velocity  profile  of  the  TF-34  turbofan 
engi ne. 

The  turbojet  test  cell,  TF-34  engine  and  augmentor  dimensions  used  in 
the  computer  model  are  shown  in  Fig.  22.  There  were  several  differences 
between  the  configuration  used  in  the  computer  model  and  the  actual  test 
cell  as  shown  in  Fig.  22: 

1.  The  computer  model  was  axi symmetric  and  therefore  required  the 
engine  to  be  located  on  the  axis  of  symmetry  of  the  test  cell. 

In  the  actual  cell  the  engine  is  approximately  0.61m  closer  to 


the  deck  than  to  the  overhead.  Thus,  it  can  be  expected  that  the 
predicted  velocity  distribution  of  the  secondary  cell  air  would 
be  somewhat  different  than  the  actual  distribution.  In  addition. 
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the  test  stand  and  supporting  deck  probably  distorted  the  velocity 
profile  at  the  augmentor  inlet. 

2.  The  exhaust  of  the  TF-34  turbofan  engine  was  modeled  as  two  con¬ 
centric  cylinders  (gas  generator  core  and  fan  exhausts).  In 
reality,  each  exit  incorporates  a  converging  nozzle.  For  simpli¬ 
city  in  the  two-dimensional  model,  the  actual  exit  diameter  of 
each  duct  was  utilized  over  the  entire  engine  length  (Figs.  23 
and  24). 

3.  In  the  computer  model  the  aft  cell  wall  was  taken  to  coincide 
with  the  augmentor  inlet.  This  wall  was  actually  approximately 
11m  farther  back.  The  effect  of  this  modification  was  to  elimi¬ 
nate  the  reci rculation  zone  above  the  augmentor  tube.  However, 
the  effects  of  this  low  velocity  reci rculation  zone  above  the 
augmentor  flow  field  have  been  shown  to  be  insignificant  [Ref.  5], 

4.  The  computer  model  utilized  a  constant  diameter  augmentor  tube 
20.5m  long  instead  of  the  varying  area  tube  construction  shown 
in  Fig.  1.  The  second  rake  was  not  positioned  farther  down¬ 
stream  to  avoid  complications  arising  from  the  divergent 
character  of  the  second  augmentor  section  and  the  flow  distur¬ 
bance  expected  to  result  from  the  water  injection  ring  and 
noise  suppression  devices.  The  augmentor  tube  length  was 
modelled  as  20.5m  versus  the  actual  5.05m  so  that  the  boundary 
conditions  at  the  exit  plane  would  not  influence  the  calcula¬ 
tions  at  the  position  of  the  second  rake. 

5.  The  test  cell  was  square  in  cross-section.  This  was  modelled  as 
a  concentric  cylinder  with  the  same  cross-sectional  area. 
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The  grid  spacing  for  the  finite  difference  numerical  solution  of  this 
problem  was  crucial.  A  30  by  30  grid  was  used  originally  without  obtain¬ 
ing  convergence.  Consequently  the  grid  was  expanded  to  a  40  by  30  to 
better  accommodate  the  much  greater  length  of  the  axial  direction.  The 


recommendation  by  Gosman 


should 


not  increase  by  more  than  a  factor  of  approximately  1.5  was  adhered  to. 

C.  DATA  INPUT 

Using  the  empirical  velocity  profiles  from  Figs.  15  through  17,  an 
average  velocity  for  each  radial  position  at  the  aft  rake  was  determined. 
These  average  velocities  and  the  centerline  velocity  were  then  plotted 
(Fig.  25)  and  graphically  integrated  to  obtain  a  one-dimensional  bulk 
velocity  at  the  aft  rake  for  each  power  setting.  An  identical  procedure 
was  followed  to  obtain  an  average  static  temperature  (Fig.  14)  at  the  aft 
rake.  These  results  were  then  used  (together  with  the  standard  properties 
of  air  and  known  augmentor  tube  cross  sectional  area)  to  calculate  the 
total  mass  flow  rate  through  the  augmentor  tube. 

Using  the  fan  and  core  RPM,  thrust,  and  engine  inlet  temperatures 
that  appear  in  Tables  II  through  IV,  fan  and  core  mass  flow  rates  for 
the  TF-34  were  obtained  from  the  General  Electric  Corp.,  Lynn,  Mass. 

It  should  be  noted  that  the  flow  rates  provided  by  General  Electric 
were  for  "an  ideal"  engine  and  may  not  exactly  apply  to  the  specific 
engine  used  in  this  investigation. 

With  the  above  flow  rates  it  was  then  possible  to  calculate  the 
mass  flow  rate  through  the  test  cell.  Then,  using  the  temperature  and 
pressure  in  the  test  cell  (and  the  known  cell  cross  sectional  area)  the 
average  one-dimensional  velocity  of  the  cell  air  was  determined. 
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In  this  investigation  no  instrumentation  was  available  for  measuring 
static  pressures  and  temperatures  at  the  TF-34  core  and  fan  tailpipe  ex¬ 
hausts.  Turbine  exit  temperature  as  a  function  of  thrust  is  presented  in 
Ref.  8. 

This  temperature  was  used  as  the  tailpipe  stagnation  temperature  in 
order  to  calculate  nozzle  exhaust  temperature.  For  the  idle  power  setting 
the  Mach  number  at  the  core  exit  was  assumed  to  be  approximately  0.5  and 
y  to  be  1.33.  Using  a  turbine  exhaust  temperature  of  1560°R  resulted  in 
a  nozzle  discharge  static  temperature  of  1498°R.  For  the  90%  and  full 
power  cases  the  nozzle  flow  was  choked  and  the  turbine  exhaust  tempera¬ 
tures  were  1694°R  and  1922°R  respectively.  This  resulted  in  exhaust 
temperatures  of  1453°R  and  1649°R  respectively.  Fan  exhaust  temperatures 
were  available  from  the  temperature  profiles  of  Ref.  8.  Actual  cross 
sectional  areas  of  both  exhaust  nozzles  were  measured  on  the  engine. 

Static  pressure  at  the  exhaust  nozzles  was  taken  to  be  the  static  pres¬ 
sure  of  the  cell  in  all  cases.  This  static  pressure  estimate  is 
accurate  for  the  fan  exhaust  in  all  cases  and  the  core  exhaust  at  the 
idle  power  setting.  However,  the  static  pressure  of  the  core  exhaust 
for  the  maximum  and  90%  power  settings  was  actually  the  critical  pres¬ 
sure. 

Using  the  above  information  the  one-dimensional  continuity  equation 
was  used  to  calculate  the  bulk  velocities  for  core  and  fan  exhausts  at 
all  power  settings.  The  computer  input  data  as  well  as  the  CPU  time 
required  to  obtain  a  converged  solution  for  each  power  setting  are  shown 
in  Table  VI. 
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POWER  SETTING 

TEMPERATURE 

(°K) 

VELOCITY 

(m/sec) 

CELL-PRESSURE 

(N/nrxl0~3) 

CPU  TIME 

(IBM  360/67) 
(min) 

CELL 

291 

0.825 

101.5 

75 

IDLE 

FAN 

344 

56.04 

CORE 

832 

59.63 

CELL 

291 

2.75 

101.37 

75 

90% 

FAN 

305 

174.2 

CORE 

807 

258.8 

CELL 

290 

3.52 

101.30 

150 

100% 

FAN 

306 

202.9 

CORE 

916 

354.1 

TABLE  VI.  COMPUTER  INPUT  DATA 
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IV.  DISCUSSION  OF  RESULTS 


A.  INTRODUCTION 

The  results  of  this  investigation  are  presented  in  the  following 
manner: 

1)  General  characteristics  and  accuracy  of  the  experimental  data 
and  sensitivity  of  the  calculated  experimental  velocity  pro¬ 
files  to  assumptions  made  for  the  unknown/unmeasured  experi¬ 
mental  parameters. 

2)  A  comparison  of  experimental  data  to  computer  model  predictions 
and  the  sensitivity  of  the  predictions  to  uncertainties  in  the 
data. 

B.  EXPERIMENTAL  DATA 

Experimental  data  are  presented  in  Tables  II  through  IV.  The  re¬ 
sulting  calculated  velocity  profiles  for  the  fore  and  aft  rakes  at  each 
power  setting  are  presented  in  Figs.  19  through  21  and  Figs.  15  through 
17,  respectively.  Figs.  10  through  12  show  the  static  pressure  profiles 
along  the  augmentor  wall  for  each  power  setting. 

The  pressure  profiles  show  only  a  small  variation  with  circumferen¬ 
tial  location.  In  all  cases  the  minimum  pressure  occurred  very  close 
to  the  inlet  plane  of  the  augmentor,  whereas  the  maximum  pressure  in 
the  instrumented  portion  of  the  augmentor  at  no  time  exceeded  atmos¬ 
pheric.  As  expected,  the  minimum  pressure  at  the  augmentor  inlet  de- 

-3  2 

creased  as  the  engine  flow  rate  increased,  going  from  -0.8x10  N/m 

-3  2 

at  the  idle  power  setting  to  -14.0x10  N/m  at  the  maximum  power  condi¬ 
tion. 
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The  velocity  profiles  at  the  aft  rake  indicate  that  the  flow  was  not 
axisymmetric  in  the  augmentor  tube.  The  most  asymmetric  condition  oc¬ 
curred  at  the  idle  power  setting.  The  idle  profile  was  also  much  more 
fully  developed  than  the  others,  indicating  that  the  flow  at  idle  was 
the  only  one  that  entered  the  diverging  section  of  the  augmentor  tube  in 
a  reasonably  well  mixed  condition.  There  also  appeared  to  be  a  general 
trend  toward  higher  mass  flow  in  the  quadrant  defined  by  the  12  and  3 
o'clock  positions. 

In  all  cases  the  velocity  profiles  at  the  forward  rake  showed  the 
expected  high  velocity  near  the  center  with  the  lower  fan  and  augmentor 
air  velocities  toward  the  wall.  It  should  also  be  noted  that  the  spread 
in  data  at  the  forward  rake  was  much  less  than  that  at  the  aft  rake. 

In  all  cases  the  engine-exit-to-augmentor-inlet  spacing  was  1.7m.  At 
the  idle  power  setting  the  augmentation  ratio  was  2.39:1.  This  ratio 
dropped  to  2.22:1  for  the  90%  power  case  and  rose  to  2.42:1  for  the  maxi¬ 
mum  power  setting. 

In  summary  it  can  be  said  that  the  static  pressure  remained  quite 
uniform  across  the  tube,  that  the  flow  at  the  aft  rake  was  not  axisymmet¬ 
ric,  and  that  at  the  front  pitot  rake  the  flow  was  quite  symmetrical. 

For  the  90%  power  setting  the  stagnation  pressure  data  at  the  aft  rake 
was  estimated  to  have  an  uncertainty  of  0.17%  which,  in  turn,  resulted  in 
a  0.25%  uncertainty  in  the  calculated  velocities.  Similarily,  the  sta¬ 
tic  pressure  readings  were  examined  and  found  to  have  an  uncertainty  of 
approximately  2%  near  the  aft  rake.  This  uncertainty  in  static  pressure 
resulted  in  approximately  a  5%  uncertainty  in  the  velocities  at  the  aft 
rake. 
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The  possibility  that  the  total  temperature  profile  at  the  aft  rake 
was  not  that  estimated  in  Figure  13  was  also  examined.  A  somewhat  dif¬ 
ferent,  but  plausible,  total  temperature  profile  was  plotted  and  is  shown 
in  Fig.  26  for  90%  power.  This  change  in  the  Ty  profile  resulted  in  a 
4.7%  increase  in  the  bulk  static  temperature  and  a  1.3%  increase  in  the 
bulk  static  temperature  and  a  1.3%  increase  in  the  bulk  velocity.  These 
changes  would  result  in  a  4.8%  decrease  in  the  test  cell  mass  flow  rate 
(assuming  engine  flow  rates  were  known). 

Finally,  the  assumption  that  the  specific  heat  ratio  was  not  1.4  at 
the  aft  rake  was  examined.  Changing  y  to  1.35  resulted  in  insignificant 
changes  in  the  cell  flow  rate. 

The  experimental  data  for  the  forward  rake  were  also  examined  to  dis¬ 
cover  how  uncertainties  in  static  pressure  and  stagnation  temperature 
profiles  would  affect  the  predicted  velocities.  Possible  errors  in  the 
extrapolated  pressure  at  the  forward  end  of  the  pitot  tubes  resulted  in 
an  uncertainty  of  only  4.5%  (Fig.  27)  in  the  velocities.  The  velocity 
profile  is  sensitive  to  the  assumed  stagnation  profile  in  much  the  same 
manner  as  discussed  above  for  the  aft  rake. 

C.  COMPUTER  MODEL  PREDICTIONS 

The  pressure  and  velocity  profiles  that  resulted  from  computer 
modelling  are  presented  with  the  actual  experimental  data  for  each 
power  setting  in  Figs.  10  through  12,  15  through  17,  and  19  through  21. 

Sensitivity  of  the  computer  model  predictions  to  specified  test 
cell  mass  flow  rate  and  engine  mass  flow  rate  were  examined.  A  T0%  de¬ 
crease  in  engine  mass  flow  rate  was  found  to  alter  the  computer  model 
results  an  insignificant  amount.  This  is  attributable  to  the  lesser 
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contribution  of  the  engine  to  the  total  flow  rate  through  the  augmentor. 

A  10%  decrease  in  the  specified  test  cell  mass  flow  rate  decreased  the 
predicted  velocity  at  the  front  rake  a  maximum  of  10%.  The  maximum  de¬ 
crease  occurred  near  the  augmentor  wall,  where  the  cell  flow  rate  had  the 
greatest  influence.  At  the  aft  rake  the  10%  decrease  in  specified  cell 
flow  rate  caused  a  maximum  decrease  in  velocity  of  7%,  again  near  the 
augmentor  wall.  This  decrease  in  specified  flow  rate  caused  the  static 
pressure  drop  to  decrease  approximately  20%.  These  results  are  shown  for 
the  90%  power  case  in  Figs.  28  and  29. 

The  theoretical  velocity  profiles  at  the  forward  rake  appear  to  be  in 
good  agreement  with  the  experimental  results.  However,  the  experimental 
results  showed  less  variations  near  the  centerline  than  predicted.  The 
predictions  underestimated  the  centerline  velocities  at  idle  and  90% 
power  and  overestimated  it  at  100%  power. 

The  experimentally  obtained  aft  rake  velocities  at  idle  power  had 
considerable  circumferential  variations.  The  theoretical  predictions, 
which  are  based  on  axisymmetric  flow,  appeared  to  be  a  reasonable  esti¬ 
mation  of  the  average  behavior.  In  both  the  90%  and  maximum  power 
cases  the  computer  model  predicted  a  good  velocity  profile  although  the 
velocities  were  considerably  greater  (up  to  63%)  than  the  experimental 
values.  This  occurred  to  a  greater  extent  at  maximum  power,  indicating 
a  possible  decrease  in  predictive  accuracy  as  mass  flow  rate  was  in¬ 
creased.  Inaccurate  specification  of  engine  flow  rates  can  significantly 
alter  the  cell  flow  rate  input  to  the  computer  model.  This  in  turn  can 
significantly  affect  the  predicted  velocity  profiles  as  discussed  above. 
Inaccurate  stagnation  profiles  can  also  contribute  to  errors  in  model 
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input  data  as  discussed  above.  However,  these  possible  inaccuracies 
cannot  account  for  the  entire  difference  in  results  obtained  at  full 


power. 


The  theoretical  pressure  profiles  (Figs.  10,  11,  and  12)  in  all 
three  cases  indicate  a  large  pressure  drop  at  the  augmentor  tube  inlet 
plane.  The  pressure  decrease  was  greater  at  the  augmentor  wall  than 
along  the  engine  centerline  .  This  theoretical  drop  in  pressure  was 
greater  than  the  experimental  pressure  drop,  the  difference  increasing 
with  an  increase  in  mass  flow  rate.  The  pressure  rise  from  the  minimum 
was  in  better  agreement  with  the  data  but  also  was  less  accurate  at 
higher  power  settings.  The  predicted  behavior  being  much  greater  than 
experimental  results  is  similar  to  the  results  of  Walters  and  Netzer 

e  and  to  the  results  of  Stevenson  and 


variable  code.  In  these  earlier 


studies  which  compared  predictions  to  experimental  data  in  a  one  eighth 
scale  test  cell,  good  agreement  in  pressure  profiles  was  obtained  for 
zero  engi ne-augmentor  spacing.  The  discrepancy  increased  with  increasing 
engine-augmentor  spacing.  One  possible  cause  of  this  disagreement  is  the 
K-e  turbulence  model.  It  is  being  used  in  this  application  for  the  shear 
layer  mixing  of  a  practically  unconfined  jet  surrounded  by  a  very  low 
velocity  flow,  both  of  which  must  then  enter  a  confined  diameter. 

Sensitivity  of  the  computer  model  predictions  to  the  magnitude  of  the 
effective  viscosity  was  examined  in  a  cursory  manner  by  both  increasing 
and  decreasing  the  calculated  values  by  a  factor  of  1.5.  The  computer 
model  predictions  were  found  to  be  fairly  insensitive  to  these  changes  in 
calculated  effective  viscosity.  For  the  case  of  a  50%  increase  in  viscosity 


38 


Figs.  30  and  31  show  less  than  5%  change  in  velocity  at  the  front  rake, 
15%  at  the  aft  rake,  and  no  greater  than  a  6%  decrease  in  static  pressure 
drop  along  the  augmentor  wall.  A  33%  decrease  in  viscosity  changed  velo¬ 
city  at  the  front  rake  less  than  5%,  10%  at  the  aft  rake  and  decreased 
the  static  pressure  drop  along  the  augmentor  wall  no  more  than  9%  (Figs. 
32  and  33).  It  appears  that  merely  changing  the  magnitude  of  the  effec¬ 
tive  viscosity  will  not  account  for  the  differences  in  the  predicted  and 
measured  pressure  profiles.  Further  investigations  of  the  two-parameter 
effective  viscosity  model  are  required. 

The  sensitivity  of  the  static  pressure  profiles  to  possible  inaccura¬ 
cies  in  model  inputs  can  account  for  much  of  the  discrepancies  at  idle 
power  but  not  for  90%  and  full  power. 

Thus,  it  appears  from  both  the  velocity  and  pressure  profile  data 
that  the  predictive  accuracy  of  the  model  decreases  with  increasing 
engine  exit  Mach  number. 

Examination  of  the  data  from  the  velocimeters  positioned  above  the 
TF-34  showed  a  test  cell  air  velocity  profile  which  increased  from 
ceiling  to  engine.  However,  it  has  been  shown  by  Speakman,  Hayes,  and 
Netzer  jRef.  5"]  that  this  condition  has  an  insignificant  effect  on  the 
predicted  velocity  and  pressure  profiles  in  the  augmentor. 

The  computer  model  appears  to  have  the  same  capabilities  and  weak¬ 
nesses  in  predicting  the  flow  in  full  scale  test  cells  as  it  did  for 
the  subscale  test  cell.  The  subscale  test  cell,  by  design,  provided 
nearly  axi symmetric  flow  in  the  augmentor  at  all  power  settings  and 
therefore  somewhat  better  agreement  was  obtained  for  predicted  velocity 
profiles  far  down  the  augmentor  tube.  Engine  flow  rates  were  also  more 
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accurately  known  in  the  subscale  tests.  The  most  severe  weakness  of  the 
model  appears  to  be  in  the  inability  to  correctly  predict  accurate  sta¬ 
tic  pressure  profiles  in  regions  where  an  expanding  high  velocity  jet  is 
captured  by  a  fixed  diameter  tube.  Velocity  profiles  and  mixing  rates 
appear  to  be  reasonably  predicted. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  two-dimensional,  primitive-variable  computer  model  proved  to  be 
quite  good  in  predicting  velocity  profiles  and  mixing  rates  in  the  aug- 
mentor  tube  of  a  full  scale  jet  engine  test  cell.  The  major  weakness  in 
this  area  was  that  the  magnitude  of  the  velocities  far  downstream  were 
somewhat  higher  than  actual  with  the  discrepancy  increasing  with  in¬ 
creasing  engine  exhaust  Mach  number  (i.e.  power  setting). 

The  greatest  drawback  to  the  computer  model  was  its  inability  to 
accurately  predict  static  pressure  profiles  in  the  augmentor  entrance 
region  where  an  expanding  high  velocity  jet  is  captured  by  a  fixed  dia¬ 
meter  tube.  This  inaccuracy  also  was  found  to  increase  with  engine 
power  level . 

If  time  and  cost  pose  less  restraints  to  future  investigators  it  is 
recommended  that  the  data  gathering  apparatus  be  expanded  to  include 
pitot-static  and  temperature  measuring  equipment  at  all  locations  where 
velocity  data  are  desired.  This  would  preclude  the  necessity  for  making 
the  assumptions  required  in  this  investigation  for  stagnation  tempera¬ 
ture  profiles. 
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Fig.  1  TF-34  Test  Configurations 


Fig.  2  Pitot  Rake  Front  View 
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Fig.  4  Pitot  Rake  Contruction 
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AUGMENTER  TUBE 


Fig.  5  Rake  Attachment 
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Fig.  6  Instrumentation 
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Fig.  7  Velocimeter  Calibration  Curves 
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Fig-  8  Pitot  Tube/Thermocouple  Labels 


Fig-  9  Static  Port  Labels 
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Fig.  10  Pressure  Profile,  Idle  Power 
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Fig.  13  Stagnation  Temperature  Profiles,  Aft  Rake 
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Fig.  14  Static  Temperature  Profiles,  Aft  Rake 
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Fig.  15  Velocity  Profile,  Aft  Rake,  Idle  Power 
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Fig.  16  Velocity  Profile,  Aft  Rake,  90%  Power 
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Fig.  17  Velocity  Profile,  Aft  Rake,  100%  Power 
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Fig.  18  Stagnation  Temperature  Profiles,  Forward  Rake 
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Fig.  19  Velocity  Profile,  Forward  Rake,  Idle  Power 
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Fig.  20  Velocity  Profile,  Forward  Rake,  90%  Power 
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Fig.  21  Velocity  Profile,  Forward  Rake,  100%  Power 
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Fig.  22  Turbojet  Test  Cell  Geometry  for  Computer  Model 


Fig.  23  TF-34  Exhaust 
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Fig.  25  Average  Velocity  Profiles,  Aft  Rake 
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Fig.  26  Modified  Stagnation  Temperature  Profile,  Aft  Rake,  90%  Power 
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Fig.  27  Sensitivity  of  Velocity  Profile  to  Uncertainties  in  Measured  Static  Pressure,  90%  Power 
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Fig.  28  Sensitivity  of  Predicted  Velocity  Profile  to  Test  Cell  Mass  Flow  Rate,  90%  Power 
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Fig.  29  Sensitivity  of  Predicted  Static  Pressure  Profile  to  Test  Cell  Mass  Flow  Rate,  90%  Power 


Fig.  30  Sensitivity  of  Predicted  Velocity  Profile  to  a  50% 
Increase  in  Effective  Viscosity,  90%  Power 
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Fig.  31  Sensitivity  of  Predicted  Pressure  Profile  to  a  50% 
Increase  in  Effective  Viscosity,  90%  Power 


Fig.  32  Sensitivity  of  Predicted  Velocity  Profile  to  a  33% 
Decrease  in  Effective  Viscosity,  90%  Power 
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Fig.  33  Sensitivity  of  Predicted  Pressure  Profile  to  a  33% 
Decrease  in  Effective  Viscosity,  90%  Power 
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